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摘 要 : 为 确定 田 块 尺度 下 探 地 雷达 对 不 同 深度 及 相 邻 反射 层 间 土壤 含水 量 的 反 演 精度 、 有 效 反 演 深度 、 
最 佳 反 演 深度 及 最 优 反 演 模型 ， 本 研究 采用 1000 MHz 中 心 频率 探 地 雷达 设备 ， 分 别 在 无 降雨 偏 干旱 土壤 和 
降雨 后 湿润 土壤 两 种 条 件 下 ， 在 选 定 农 田 区 域 基于 共 中 心 点 法 采集 雷达 波 数据 ， 提 取 有 效 地 表 波 与 反射 波 
数据 ， 通 过 双 曲 线 拟 合法 分 别 获 取 不 同 深度 反射 层 雷达 波 的 传播 速度 ， 计 算得 出 土壤 的 相对 介 电 常数 ， 最 
后 根据 土壤 体积 含水 量 和 介 电 常数 之 间 的 经 验 模 型 计算 获得 不 同 深度 的 土壤 体积 含水 量 。 通 过 Topp Roth, 
Herkelrath 和 Ferre 四 种 经 验 模型 分 别 进行 土壤 体积 含水 量 反 演 测定 ， 同 时 以 利用 烘 干 法 获取 的 代表 性 测定 
土壤 含水 量 实测 值 为 指标 进行 精度 验证 。 田 间 试 验 结果 表明 ，1000 MHz 探 地 雷达 的 有 效 反 演 深度 范围 为 
0—50 em; 土壤 偏 干旱 和 偏 湿 润 条 件 的 最 佳 反 演 深度 分 别 为 50 cm 和 40 cm; Roth 模 型 相关 性 最 好 ， 决 定 系 
数 尼 最 高 为 0.750， 且 Roth 模 型 反 演 土壤 含水 量 值 最 稳定 ， 在 土壤 偏 干旱 和 偏 湿 润 条 件 下 均 方 根 误差 (Root 
Mean Square Error, RMSE) 平均 值 分 别 为 0.0401 mzym 和 0.0335 m*/m*; 相对 误差 (Relative Error, RE) 最 低 
为 3.0%。 探 地 雷达 具备 对 定量 深度 田间 土壤 含水 量 快速 、 精 准 监测 的 能 力 ， 但 其 反 演 模型 需 根据 不 同 土壤 
类 型 等 条 件 进行 相应 参数 校正 。 
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1 引言 量 平衡 和 降水 分 布 具 有 十 分 重要 的 意义 ""。 土 壤 
中 水 分 含量 的 变化 会 对 作物 生长 状况 产生 直接 影 

从 全 球 水 文 循 环 来 看 ， 土 壤 水 分 总 量 较 小 是 ” 响 。 特 别 在 田 块 尺度 上 ， 土 壤 含水 量 是 作物 产量 

仅 占 全 部 水 分 总 量 的 0.05% 左 右 ， 但 其 对 全 球 能 “和 品质 的 主要 控制 因素 中 。 因 而 ， 有 效 精准 地 测 
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定 田 间 土 壤 水 分 含量 对 作物 长 势 监 测 和 产量 预测 
具有 重要 意义 。 

现 阶段 ， 研 究 人 员 基 于 不 同 的 需要 ， 对 于 土 
培 含 水 量 的 测定 已 有 较 多 的 技术 方案 。 传 统 测 定 
土壤 含水 量 的 方法 主要 有 烘 干 法 、 中 子 法 以 及 时 
域 反射 法 (Time Domain Reflectometry, TDR) 
等 。 这 些 方法 能 较 准 确 地 测量 指定 位 置 的 土壤 含 
水 量 ， 但 在 实验 操作 中 对 田 十 具有 破坏 性 。 卫 星 
遥感 技术 的 出 现 打 破 了 区 域 尺度 大 面积 土壤 水 分 
动态 监测 的 壁垒 ”， 人 研究 人 员 通 过 建立 经 验 或 半 
经 验 模 型 来 提高 土壤 含水 量 的 反 演 精度 ， 但 卫星 
遥感 易 受 植被 覆盖 、 天 气 条 件 等 影响 “， 并 且 和 针 
对 田 块 尺度 下 不 同 深 度 土 壤 含水 量 的 监测 ， 卫 星 
遥感 数据 的 分 辨 率 、 保 障 率 和 精度 有 待 提 高 ”。 
因此 ， 田 块 尺度 定量 深度 土壤 含水 量 的 高 精度 、 
快速 无 损 测 量 成 为 目前 研究 的 难点 。 

探 地 雷达 (Ground Penetrating Radar, GPR) 
是 近年 来 新 兴 的 一 种 土壤 水 分 监测 技术 手段 ， 与 
众多 十 壤 含水 量 探测 方法 相 比 ， 探 地 雷达 具有 操 
作 简 单 、 分 辩 率 高 、 探 测 深度 大 、 可 重复 性 以 及 
对 地 下 非 侵入 性 测量 等 优点 。 由 于 控 地 雷达 具有 
上 述 优势 量 在 测量 过 程 中 无 需 扰动 作物 与 土壤 ， 
所 以 在 田 块 尺度 深层 土壤 含水 量 的 测量 方面 表现 
出 很 大 的 潜力 。 综 合 国内 外 探 地 雷达 研究 现状 ， 
目前 基于 探 地 雷达 和 雷达 波 信号 分 析 的 土壤 水 分 
快速 无 损 监测 研究 主要 集中 在 以 下 3 个 方面 “*。 

(1) 快速 测定 土壤 表层 含水 量 。Grote 等 ” 
使 用 900 MHz 和 450 MHz 探 地 雷达 测量 美国 加 州 
葡萄 园 表 层 土壤 含水 量 ， 雷 达 反 演 值 与 烘 干 法 实 
测 值 相关 性 较 好 ,证 明了 探 地 雷达 地 表 波 可 以 快 
XE. 、 无 破坏 性 地 提供 大 面积 土壤 表层 含水 量 。 
Galagedara 等 " 使 用 探 地 雷达 和 6 个 不 同 长 度 的 
时 域 反 射 仪 ， 在 沙壤土 上 设计 排水 试验 探究 探 地 
雷达 的 有 效 反 演 深 度 ， 结 果 表 明 100 MHz 探 地 雷 
达 在 0~-0.1 m 深 度 反 演 结 果 较 精确 。 

(2) 不 同 地 质 类 型 的 土壤 含水 量 反 演 。Al- 
geo 等 … 使 用 探 地 雷达 振幅 分 析 法 监测 粘性 土壤 
中 土壤 含水 量 对 灌溉 的 响应 变化 ， 试 验 结果 表明 


该 方法 可 以 获取 富 含 粘 粒 土 壤 中 地 下 含水 率 的 空 
间 分 布 信息 。 王 前 锋 等 ”利用 宽 角 反 射 与 折射 
法 (Wide Angle Reflection and Refraction, 
WARR) 和 固定 间距 法 (Fixed Offset, FO) 相 结 
合 证 明了 200 MHz 中 心 频率 天 线 在 沙漠 地 区 地 波 
法 的 最 佳 天 线 间 距 为 1.9 m， 有 效 深度 为 20 cm, 
偏差 仅 为 0.015 mYm’。 曹 棋 等 ' ”利用 探 地 雷达 
地 波 法 证 明了 60 MHz 天线 在 红壤 地 区 反 演 0— 
40 cm 土 层 的 土壤 体积 含水 量 时 精度 最 高 Roth 
经 验 模型 更 适合 于 红壤 地 区 ， 均 方 根 误差 (Root 
Mean Square Error, RMSE) 为 0.022 m"/m', 

(3) 土壤 含水 量 反 演 模 型 建立 及 参数 校正 。 
李 草 君 和 钟 若 飞 "利用 探 地 雷达 探究 土壤 介 电 
常数 与 雷达 回 波 振幅 之 间 的 变化 规律 ， 建 立 了 探 
地 雷达 波 振幅 与 土壤 介 电 常 数 之 间 的 定量 关系 ， 
模型 反 演 精度 较 高 ，RMSE 为 0.014 m/m. Cao 
A US 利用 探 地 雷达 探究 了 地 表 至 1 m 深 的 三 维 
土壤 水 分 变化 以 及 降雨 事件 对 土壤 水 分 空间 动态 
的 影响 ， 并 每 间隔 20 cm 建立 了 介 电 常数 与 土壤 
体积 含水 量 的 反 演 模型 。 

目前 大 部 分 研究 主要 侧重 在 不 同 地 质 类 型 的 
土壤 含水 量 反 演 以 及 经 验 模型 的 校正 ， 基 于 田 块 
尺度 的 不 同 深 度 土壤 层 含水 量 反 演 的 相关 研究 较 
少 ， 实 验 数 据 尚 不 充分 。 然 而 ， 在 田 块 尺度 下 ， 
不 同 深度 土壤 层 水 分 对 作物 的 生长 会 产生 相应 不 
同 的 影响 "” ， 能 否 快 速 、 无 损 地 对 不 同 深度 土 
培 层 含水 量 进行 监测 对 于 田间 精准 灌溉 、 作 物 生 
长 监测 、 作 物品 质 与 产量 预测 等 方面 起 到 重要 
作用 。 

因此 ， 为 确定 田 块 尺度 下 探 地 雷达 对 不 同 深 
度 及 相 邻 反射 层 间 土壤 含水 量 的 反 演 精度 、 有 效 
反 演 深度 、 最 佳 反 演 深 度 及 最 优 反 演 模 型 ， 本 研 
究 采 用 1000 MHz 中 心 频率 探 地 雷达 设备 ， 分 别 
在 无 降雨 偏 干 旱 土壤 和 降雨 后 湿润 土壤 条 件 下 的 
标准 农田 环境 ， 反 演 不 同 深度 土壤 体积 售 水 量 ， 
并 与 对 应 取样 土壤 烘 干 法 实测 值 进行 对 比 验 证 。 
同时 ， 利 用 雷达 波 从 地 表 传 播 至 不 同 深度 反射 层 
的 平均 速度 和 迪克 斯 (Dix) 公式 ,探究 探 地 雷 
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达 对 相 邻 反射 层 间 土壤 水 分 的 反 演 精度 。 
2 材料 与 方法 


2.1 研究 区 概况 和 试验 设备 


研究 区 位 于 北京 市 顺义 区 大 孙 各 庄 镇 中 国 
农业 科学 院 农 业 环境 与 可 持续 发 展 研究 所 试验 
基地 。 顺 义 区 年 均 降 雨量 约 625 mm， 全 年 降水 
的 75% 集 中 在 夏季 。 试 验 选取 人 研究 区 的 50 mX 
20 m 试 验 农田 ， 农 田 土壤 以 潮 土 、 褐 土 为 主 ， 
作物 常年 为 小 麦 、 玉 米 轮 作 。 试 验 期 间 农 田 无 作 
物种 植 ， 雷 达 移 动 路 径 为 裸露 地 表 且 地 表 较 平 
整 ， 适 宜 进行 雷达 试验 等 。 

试验 采用 加 拿 大 Sensor & Software 公 司 生 产 
的 Pulse EKKO PRO 型 探 地 雷达 ， 由 主机 、 发 射 
天 线 (Transmitter，T) 、 接 收 天 线 (Receiver, 
R) 、 电 池 及 其 连接 装置 构成 ， 其 中 两 个 天 线 频 
率 均 为 1000 MHz。 其 它 试 验 设备 为 50 m 卷 尺 以 
及 烘 干 法 采集 土 样 所 需 的 土 钻 、 环 刀 、 标 签 、 铝 
盒 等 。 

2.2 试验 方法 
2.2.1 探 地 雷达 测量 土壤 含水 量 的 原理 

探 地 雷达 发 射 天 线 向 外 发 射 高 频 电磁 波 ， 部 
分 电磁 波 通 过 空气 直接 从 发 射 天 线 传播 到 接收 天 
线 ， 这 些 电 磁 波 称 为 空气 波 (Air Wave, AW). 
地 表 波 (Ground Wave, GW) 沿 着 土壤 表面 传 
播 到 接收 天 线 ， 甚 路程 可 近似 看 作 两 个 天 线 之 间 
的 间距 。 

反射 波 (Reflected Wave, RW) 或 折射 波 穿 
过 地 表 到 达 地 下 介质 ， 当 电磁 波 遇 到 电 性 差异 大 
的 分 界面 、 目 标 体 或 者 具有 不 同 介 电 常数 的 材料 
时 ， 探 地 雷达 电磁 波 的 振幅 、 频 率 和 相位 都 会 发 
生变 化 ， 从 而 发 生 散射 和 反射 ， 并 从 不 同 的 界面 
返回 接收 天 线 。 随 后 ， 接 收 天 线 接收 各 种 电磁 
波 ， 根 据 探 地 雷达 电磁 波 的 双 程 走时 、 相 位 、 振 
幅 等 参数 ， 可 以 确定 地 层 结构 或 目标 体 的 埋 深 、 
位 置 等 "7 。 

电磁 波 在 介质 中 传播 速度 取决 于 光 在 真空 中 


传播 速度 、 相 对 介 电 常数 (e) 和 相对 磁 导 
XC. 雷达 波 在 非 饱 和 带 土 壤 中 的 传播 速度 与 
土壤 的 介 电 常 数 密切 相关 。 常 温 下 ， 土 壤 固 相 的 
相对 介 电 常数 大 约 是 3~10， 空 气 的 相对 介 电 党 
数 是 1， 淡 水 的 相对 介 电 常数 是 81， 淡 水 的 相对 
介 电 常数 远大 于 空气 和 干 土 壤 的 相对 介 电 常数 。 
因此 ， 非 饱和 土壤 的 介 电 稼 数 主要 依赖 于 土壤 含 
水 量 “”。 根据 雷达 波 传播 速度 与 介 电 常数 之 间 
的 关系 式 (公式 (1) )， 可 以 计算 出 土壤 的 相对 
介 电 常 数 e : 

e-(£] (1) 


Vy 

其 中 ,vy 为 雷达 波 在 土壤 中 的 平均 速度 ，m/ns; 
c 为 自由 空间 中 电磁 波 传播 速度 ，0.3 m/ns。 

因此 ， 可 以 利用 共 中 心 点 法 雷达 剖面 图 像 提 
取出 雷达 波 在 土壤 中 的 传播 速度 ， 通 过 公式 (1) 
获得 土壤 的 相对 介 电 常数 ， 最 后 根据 土壤 体积 含 
水 量 和 介 电 常数 之 间 的 经 验 模 型 或 半 理 论 关系 式 
获得 土壤 体积 含水 量 。 

为 描述 介 电 常数 与 土壤 体积 含水 量 0 之 间 
的 关系 ， 研 究 者 们 已 经 提出 了 很 多 经 验 模型 ， 其 
中 Topp 模型 2 Roth {R AI °° 、Herkelrath Bi 
型 3 pL Ferre 8 ^9 因 其 变量 少 、 表 达 式 简 
洁 、 相 关 性 好 而 被 广泛 应 用 ， 故 本 研究 采用 这 4 
个 模型 对 土壤 含水 量 进行 反 演 。 模 型 表达 式 见 
公式 (2) ~ 公式 (5)。 


Topp 模型 : 

= -0.053 + 0.0292¢ - 0.00055¢? + 
0.00000432? (2) 

Roth 模型 ; 

0 = -0.0728 + 0.04482 - 0.001952? + 
0.00003612? (3) 

Herkelrath 模型 . 

6 = -0.1626+ 0.1117 Ve (4) 

Ferre 模 型 ; 

6 = -0.1841 + 0.1181Ve (5) 


2.2.2 Hab AE 
共 中 心 点 法 (Common Midpoint, CMP) 是 
在 待 测 点 确定 后 ， 将 探 地 雷达 发 射 天 线 与 接收 天 
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线 以 待 测 点 为 中 心 位 置 水 平 放置 ， 并 同时 向 相反 
的 方向 等 距离 移动 的 测量 方法 (图 1 a) )。 共 
中 心 点 法 可 以 有 效 地 提高 探 地 雷达 图 像 效 果 。 当 


共 中 心 点 法 (CMP) 


(a) 共 中 心 点 法 示意 图 


反射 层 


(Cc) 反射 疲 双 程 旅行 轨迹 


使 用 共 中 心 点 法 收集 雷达 数据 时 ， 由 于 中 心 位 置 
不 变 ， 雷 达 波 的 登 加 提高 了 地 下 雷达 反射 的 信 噪 
比 ， 从 而 使 雷达 剖面 图 像 更 清晰 Is 
空气 波 AW 


地 表 波 GW 
反射 波 RW 


雷达 波 双 程 旅行 时 间 /ns 


0 02 04 0.6 0.8 1.0 
天 线 间 距 /m 


(b) 共 中 心 点 法 雷达 图 像 


图 1 共 中 心 点 法 示意 图 和 雷达 波 速度 图 


Fig. 1 Schematic diagram of common midpoint(CMP ) method and radar-wave velocity diagram 


试验 中 雷达 收发 天 线 每 次 分 别 向 左右 移动 
0.01 m， 移 动 距离 的 总 长 度 为 2m。CMP 雷 达 图 
像 是 由 100 道 雷达 数据 组 成 (图 1 (b) )， 每 一 
道 雷达 数据 表示 在 相同 的 天 线 间 距 下 ， 雷 达 波 在 
介 电 常数 不 同 的 反射 层 发 生 反 射 ， 分别 以 不 同 的 
双 程 旅行 时 间 被 接收 天 线 获 取 ， 并 在 主机 中 以 图 
像 形式 记 录 。100 道 雷达 数据 组 成 雷达 图 像 中 的 
若干 条 灰色 曲线 ， 从 图 像 中 可 以 清晰 拾取 地 表 波 
与 空气 波 。 由 于 表层 土壤 反射 层 较 浅 ， 且 雷达 发 
射 天 线 与 接收 天 线 距离 较 小 ， 因 此 地 表 波 在 雷达 
图 像 中 表示 为 一 条 直线 ， 其 斜率 代表 地 表 波 的 波 
速 。 空 气 波 是 雷达 波 从 发 射 天 线 经 由 空气 直接 由 
接收 天 线 接收 ， 由 于 电磁 波 在 空气 中 的 速度 大 于 
在 土壤 中 的 速度 ， 因 而 最 先 被 接收 天 线 获取 ， 所 
以 在 雷达 图 像 中 以 第 一 条 直线 表示 ， 理 论 上 其 速 
度 值 接近 0.3 m/ns。 本 研究 雷达 图 像 未 能 完全 将 
空气 波 捕 获 ， 文 中 的 空气 波 仅 代表 雷达 波 通 过 与 
空气 接触 的 浅 表层 土壤 的 雷达 图 像 ， 但 本 人 研究 并 
未 应 用 空气 波 的 相关 性 质 ， 因 此 对 试验 结果 不 造 


影响 。 

雷达 图 像 中 地 表 波 以 下 的 每 条 曲线 代表 介 电 
稼 数 不 同 的 层 间 反 射 ， 即 为 反射 波 ， 其 反射 曲线 
的 曲率 代表 反射 波 的 传播 速度 。 由 于 土壤 反射 层 
相 比 地 表 深 度 较 大 ， 因 而 不 能 忽略 反射 层 深度 对 
雷达 波 传播 距离 的 影响 ， 即 不 能 将 雷达 波 运动 过 
程 视 为 雷达 天 线 间距 。 图 1 (c) 为 反射 波 在 土壤 
层 中 的 传播 过 程 ， 其 中 4 表示 反射 层 距 离 地 表 的 
EARE, m; x 代表 雷达 天 线 间 的 水 平 距 离 ， 
m; y 代 表 反 射 波 在 土壤 中 传播 的 平均 速度 ， 
m/ns; ! 代 表 反 射 波 双 程 旅行 时 间 ，ns。 根 据 
图 1 (c) 雷达 波 双 程 旅行 时 间 与 天 线 间 的 水 平 距 
离 的 关系 可 知 : 


(L«] -(Ls] oe © 


2 2 
整理 得 : 

2 2 

E recs cd (7) 
E m 

y 
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又 因为 x>0， 所 以 雷达 反射 波 的 图 像 在 xt 坐 
标 系 下 是 焦点 在 t 轴 的 双 曲 线 的 一 半 ( 取 1 轴 正 半 
轴 部 分 )。 因 此 ， 本 研究 在 获取 测 点 不 同 深度 反 
射 层 的 反射 波 速度 时 ， 应 用 双 曲 线 拟 合法 直接 从 
软件 读 取 ， 而 不 通过 计算 反射 波 图 像 近似 直线 的 
斜率 获得 ， 双 曲线 拟 合 图 像 如 图 1 (D) 反射 波 
所 示 。 
通过 上 述 方法 获取 的 速度 值 代表 雷达 波 从 地 
表 传播 至 不 同 深度 反射 层 的 平均 速度 值 ， 其 反应 
了 地 表 至 不 同 深 度 反 射 层 的 平均 含水 量 。 若 想 分 
别 反 演 不 同 深度 相 邻 反射 层 间 土 壤 体积 含水 量 ， 
需 获取 雷达 波 在 土壤 相 邻 反射 层 中 移动 的 速度 ， 
因此 本 研究 尝试 采用 Dix 公 式 ”将 雷达 波 至 反 
射 层 的 平均 速度 转换 为 雷达 波 在 相 邻 反射 层 间 移 
动 的 层 速度 ， 其 表达 式 如 公式 (8). 


2 


2 2 
Li soil, n trw, n v soil, n — ifaw, n-1 
Vine, n ( 8 ) 


lawn 7 Érw,n-1 


其 中 ， von JP n RZ EBOEIE 
XE, m/ns; ter, NE n ATE K BC BOE T EST 
B], ns; v1 为 第 n-1 个 反射 层 之 上 的 平均 波 
速 ，m/ns; trw ,1 为 第 n-1 个 反射 层 的 反射 波 旅 
行 时 间 ，ns。 


2.3 数据 采集 与 处 理 


试验 数据 采集 共 分 为 雷达 数据 和 不 同 深度 土 
壤 样 本 两 部 分 。 其 中 ， 雷达 数据 和 土壤 样本 采集 
均 两 次 ， 第 一 次 数据 获取 在 2021 年 6 月 28 日 至 7 
月 1 日 ， 试验 期 间 无 降雨 事件 ， 最 高 温度 为 
33 °C; 第 二 次 数据 获取 在 2021 年 10 月 12 日 至 


烘 烤 12 小 时 ， 待 烘箱 内 冷却 至 室温 后 取出 ， 称 
重 并 计算 不 同 次 度 土壤 层 体积 含水 量 。 由 于 土 外 
取样 过 程 是 在 同一 测 点 的 连续 10 cm 深度 取样 ， 
此 分 别 将 同一 测 点 不 同 深 度 的 土壤 层 体积 
含水 量 取 平 均值 计算 0—10 cm, 0~20 em, 
0~30 cm, 0~40 cm, 0~50 cm 和 0~-60 cm 土壤 
体积 含水 量 。 

探 地 雷达 CMP 数据 采 集 以 烘 干 法 取样 点 作 
为 中 心 点 ， 无 降雨 偏 干旱 土壤 和 降雨 后 湿润 土壤 
条 件 下 各 采集 15 幅 雷达 图 像 。 使 用 探 地 雷达 数 
据 处 理 软件 EKKO_ Project 对 雷达 数据 进行 增益 、 
滤波 等 处 理 ， 有 利于 雷达 波 的 提取 和 识别 。 地 表 
波 的 速度 由 地 表 波 所 在 直线 的 斜率 获取 ， 反 射 波 
的 速度 由 双 曲 线 拟 合法 从 软件 中 直接 读 取 ， 将 获 
取 的 速度 值 分 别 带 入 公式 〈1)， 进 而 求 得 地 表 至 
AN FIRE EERTE ER, BA e HI 
模型 (公式 (2) ~ 公式 (5) ) 可 以 获得 相应 的 
土壤 体积 含水 量 模型 预测 值 。 将 4 个 模型 的 不 同 
深度 土壤 体积 含水 量 模型 预测 值 分 别 与 对 应 深度 
土壤 烘 干 法 实测 值 进行 精度 评价 ， 依 据 精 度 评价 
结果 判断 探 地 雷达 对 不 同 深度 土壤 含水 量 反 演 精 
度 、 有 效 反 演 深度 、 最 佳 反 演 深度 以 及 1000 MHz 
探 地 雷达 田 块 尺度 最 优 反 演 模型 。 采 用 RMSE、 
R^. #HXtiRZ (Relative Error, RE) 和 绝对 误差 
(Absolute Error, AE) 对 试验 结果 进行 精度 评 
价 。 同 时 ， 结合 已 从 雷达 图 像 中 获取 的 地 表 波 和 
反射 波 的 速度 值 ， 利 用 Dix 公 式 将 雷达 波 的 平均 
速度 转换 为 层 速 度 ， 进 而 结合 最 优 反 演 模 型 获取 


15 日， 试验 前 两 天 内 有 少量 降雨 ， 田 间 土 壤 湿 
润 。 并 且 两 次 取样 位 置 相同 ， 雷 达 数 据 与 土壤 样 
本 数据 一 一 对 应 。 

在 试验 农田 间隔 5 m 平 行 布置 三 条 50 m 长 的 
测 线 ， 由 测 线 端点 开始 每 间隔 8m 插 上 地 标 ， 作 
为 土壤 烘 干 法 取样 点 。 用 环 刀 与 土 钻 在 每 个 取样 
点 处 取 0~10 cm, 10~20 cm, 20~30 cm, 30~ 
40 cm, 40~50 cm Fil 50~60 cm 六 个 不 同 深度 土 
层 的 样本 ， 密 封 并 做 好 标记 后 放置 烘箱 中 105 °C 


不 同 深度 相 邻 反射 层 间 的 土壤 体积 含水 量 模型 反 
演 值 。 通 过 将 模型 反 演 值 与 土壤 含水 量 烘 干 法 实 
测 值 进行 精度 评价 ， 探 究 探 地 雷达 对 不 同 深 度 相 
邻 反射 层 间 土壤 含水 量 的 反 演 精度 。 


RMSE = [= ce ~ Bond) (9) 


R=1- > Tun 区 Eu (10) 


> 
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3 结果 与 分 析 


3.1 土壤 体积 含水 量 及 雷达 有 效 反 演 深度 


两 次 试验 CMP 测 点 的 0~60 cm 不 同 深度 土 
壤 体 积 含水 量 烘 干 法 实测 值 如 图 2 所 示 。 图 像 的 
中 心 点 表示 相同 深度 土壤 样本 体积 含水 量 实测 值 
的 平均 值 ， 中 心 点 处 有 色 部 分 的 垂直 长 度 表 示 所 
在 深度 范围 土壤 体积 舍 水 量 的 标准 差 。 图 2 (a) 
中 不 同 深 度 土壤 体积 含水 量 的 标准 差 依 次 为 
0.03、0.02、0.02、0.01、0.01、0.01; 图 2 (b) 


o 
[m 
E 


> 
m 
N 


o 
b 
© 


ERER EIK ENm m’) 


o 
iS) 
oo 


024 1 : H A : : 
0-10 0-20 0-30 0-40 0-50 0-60 
深度 /cm 


(a)2021 年 6 月 28 日 至 7 月 1 日 无 降雨 偏 干旱 条 件 土 壤 体 积 含水 量 


实测 值 


中 标准 差 依 次 为 0.03、0.02、0.03、0.02、 
0.02, 0.02. 

从 图 中 可 以 看 出 ， 在 无 降雨 条 件 下 ， 田 间 表 
层 0~10 cm 土壤 含水 量 偏 低 ， 样 本 最 低 土 壤 含 水 
量 为 0.24。 随 着 采样 深度 逐渐 增加 ， 土 壤 合 水 量 
逐渐 升 高 ， 在 60 cm 处 最 大 值 为 0.31。 本 次 试验 
前 的 降雨 主要 对 表层 土壤 水 分 产生 影响 ，0~ 
10 cm 深度 范围 土壤 含水 量 最 大 值 为 0.36。 由 于 
环 刀 采样 时 表层 土壤 水 分 没有 完全 下 渗 ， 因 而 随 
着 深度 的 增加 ， 土 壤 含 水 量 整体 呈 降 低 趋 势 。 


0~10 0~20 0~30 0~40 0~50 0~60 


深度 /cm 


(b)2021 年 10 月 12 日 至 15 日 降雨 后 偏 湿 润 条 件 土壤 体积 含水 量 
实测 值 


图 2 无 降雨 偏 干 旱 和 降雨 后 偏 湿 润 土壤 条 件 下 CMP 测 点 0~-60 cm 不 同 深度 土壤 体积 含水 量 实测 值 


Fig.2 Soil volumetric water contents from 0 to 60 cm of CMP test points under arid soil and moist soil conditions 


图 3 是 两 次 雷达 数据 采集 的 CMP 法 雷达 图 
像 。 其 中 图 3 (a) 一 图 3 (c) 为 无 降雨 条 件 的 
CMP AIR; 图 3 (d) ~ 图 3 (f) 为 降雨 事件 后 
的 CMP 图 像 。 为 更 清楚 地 获取 不 同 采 样 条 件 下 
CMP 图 像 特点 ， 图 3 所 选 三 个 测 点 均 不 在 测 线 同 
一 行 、 同 一 列 ， 且 两 次 试验 CMP 和 采集 位 置 均 相 
同 。 从 图 像 中 可 以 很 容易 发 现 ， 两 次 雷达 图 像 的 
最 大 时 延 相差 较 小 ， 表 明 雷 达 波 在 两 次 试验 中 最 
大 探测 深度 非常 接近 。 但 从 雷达 波 获 取 情 况 来 
看 ， 在 无 降雨 条 件 下 ， 雷 达 波 通过 不 同 深度 反射 
层 所 产生 的 反射 波 分 层 更 明显 。 而 降雨 后 受 土 壤 
中 下 渗水 分 影响 ,雷达 图 像 在 相 邻 反射 层 间 容 易 
出 现 杂 波 ， 具 体 表现 为 短小 曲线 ， 随 着 天 线 间距 


的 增加 ， 其 杂 波 也 随 之 消失 。 

图 4 是 两 次 试验 利用 CMP 法 获取 的 测 点 处 
0 一 50 cm 土壤 相对 介 电 常数 分 布 图 。 中 心 点 表示 
相同 深度 范围 的 土壤 相对 介 电 常数 平均 值 ， 中 心 
点 处 有 色 部 分 垂直 长 度 表 示 相 对 介 电 常数 的 标准 
差 。 由 于 30 个 测 点 的 CMP 图像 几乎 在 60 cm 
处 无 法 获取 明显 反射 层 ， 并且 雷达 图 像 在 SO~ 
60 cm 深度 范围 没有 明显 呈 双 曲线 状 的 反射 波 ， 
因此 本 研究 后 续 土 壤 含 水 量 反 演 精度 仅 讨论 至 
50 cm, 

图 中 0~10 cm 深度 范围 的 土壤 介 电 常数 由 地 
表 波 斜率 求 得 ，0~-20 cm 至 0~50 cm 土壤 介 电 党 
数 由 反射 波 双 曲线 拟 合法 在 软件 中 直接 读 取 。 两 
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a)2021 年 6 月 28 日 至 7 月 1 日 无 降雨 偏 干旱 条 件 土壤 相对 介 电 常数 (b)2021 年 10 月 12 日 至 15 日 降雨 后 偏 湿润 条 件 土 壤 相 对 介 电 常数 
图 4 无 降雨 偏 干旱 和 降雨 后 偏 湿 润 土壤 条 件 下 CMP 测 点 0~50 cm 土壤 相对 介 电 常数 值 


Fig. 4 Soil relative dielectric constant from 0 to 50 cm of CMP test points under arid soil and moist soil conditions 
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种 土壤 水 分 条 件 下 ， 不 同 深度 土壤 的 介 电 常数 值 
在 分 布 图 中 分 层 明显 ， 无 降雨 条 件 的 土壤 介 电 常 
数 与 土壤 含水 量 实测 值 变化 趋势 一 臻 ，0~10 cm 
土壤 介 电 常数 偏 低 ， 随 着 深度 的 增加 介 电 常数 
值 也 逐渐 增 大 。 与 无 降雨 条 件 相 比 ， 降 雨 后 0— 
20 cm 深度 范围 的 土壤 介 电 常数 明显 偏 高 ， 但 与 
降雨 后 土壤 含水 量 实测 值 变化 趋势 不 同 的 是 ， 
0—20 cm £ 0—50 cm 土壤 相对 介 电 常数 随 深度 的 
增加 逐渐 增 大 。 
相同 深度 范围 土壤 的 介 电 常数 在 +3 小 范围 
内 波动 ， 说 明 地 下 相同 深度 的 土壤 反射 层 分 布 不 
均匀 。 因 此 在 使 用 双 曲 线 拟 合法 处 理 雷 达 数 据 
时 ， 其 反射 层 深度 值 变化 幅度 在 2 cm 左右 均 视 
为 正常 值 ， 且 对 于 在 某 一 深度 范围 CMP 图 像 中 
未 有 明显 反射 波 的 测 点 数据 ， 在 进行 相关 性 分 析 
VA Be Jit BE RMSE 等 精度 验证 时 为 了 对 整体 精度 验 
证 结果 不 造成 较 大 影响 ， 在 计算 时 应 将 图 像 中 无 
反射 波 的 深度 范围 对 应 土壤 水 分 实测 值 同 时 
删除 。 
通过 将 不 同 深度 土壤 体积 含水 量 实测 值 与 对 
应 深度 范围 土壤 相对 介 电 常数 进行 相关 性 分 析 发 
现 ( 表 1)。 两 次 试验 0~10 cm 深度 范围 的 相关 
系数 均 较 低 ， 表 明 两 次 试验 中 雷达 地 表 波 对 表层 
土壤 介 电 常 数 的 反 演 精度 偏 低 ， 其 原因 可 能 是 由 
于 1000 MHz 探 地 雷达 的 地 表 波 有 效 代 表 深 度 的 
不 确定 性 ， 导 致 实际 土壤 样本 采样 深度 与 地 表 波 
代表 深度 不 匹配 ， 从 而 造成 误差 。 将 两 次 不 同 土 
壤 水 分 条 件 下 的 相关 系数 对 比分 析 发 现 ， 降 雨 对 
相关 系数 的 影响 较 小 ， 相 关系 数 均 随 着 土壤 深度 
的 逐渐 增加 而 逐渐 增 大 ， 并 在 0~40 cm 处 达到 最 
大 值 。 同 时 ， 两 次 试验 结果 在 0~50 cm 深度 范围 
均 有 较 好 的 相关 性 ， 因 而 相关 性 分 析 结 果 表 明 
1000 MHz 探 地 雷达 在 本 次 试验 田 块 的 有 效 反 演 
深度 为 50 cm. 


3.2 田 块 尺度 最 优 模型 及 最 佳 反 演 深度 


通过 雷达 波 图 像 获得 的 土壤 相对 介 电 常数 结 
合 Topp 模 型 、Roth 模 型 、Herkelrath 模 型 和 Ferre 


表 1 无 降雨 偏 干 旱 和 降雨 后 偏 湿润 土壤 条 件 下 不 同 深度 
范围 土壤 介 电 常数 与 烘 干 法 实测 体积 含水 量 值 之 间 的 相 
关 性 分 析 (n = 30) 


Table 1 The correlation analysis between soil relative dielec- 
tric constant and volumetric water content under arid soil 


and moist soil conditions in different depth ranges(n = 30) 


u 相关 系数 (r) 
土壤 水 分 条 件 
0—10 cm 0~20 cm 0—30 cm 0~40 cm 0—50 cm 
偏 干旱 0.554 0.593 0.707 0.854 0.740 
偏 湿润 0.475 0.686 0.782 0.847 0.750 


模型 ， 即 可 得 到 土壤 体积 含水 量 的 模型 预测 值 。 
图 $ 是 1000 MHz 探 地 雷达 分 别 在 两 种 土壤 水 分 
条 件 下 ， 不 同 深 度 土壤 烘 干 法 实测 值 与 4 个 模型 
反 演 值 的 相关 性 对 比 图 。 其 中 ， 图 5 (a) 代表 无 
降雨 条 件 下 的 雷达 反 演 精度 ， 图 中 灰色 十 字 代 表 
在 有 效 反 演 深度 范围 4 个 模型 的 最 小 RMSE 平 均 
值 ， 其 所 在 位 置 代表 在 该 深度 范围 具有 最 小 
RMSE 的 模型 。 

在 无 降雨 条 件 下 ，Topp 模 型 的 相关 性 随 着 


: Herkelrath 


5 Herkelrath 

oll 一 一 Ferre . ol 一 一 Ferre : 

“0-10 0-20 0-30 0-40 0-50 “0-10 0-20 0-30 0-40 0-50 
深度 /cm 深度 /cm 


(a) 无 降雨 土壤 偏 干 旱 条 件 土壤 (b) 降 雨 后 土 壤 偏 湿润 条 件 土 
含水 量 模型 反 演 值 与 实测 值 相 ” 壤 含水 量 模 型 反 演 值 与 实测 值 
关 性 分 析 相关 性 分 析 

图 5 无 降雨 偏 干旱 和 降雨 后 偏 湿润 土壤 条 件 下 0~50 cm 5 

个 深度 范围 4 种 模型 预测 值 与 烘 干 法 土壤 体积 含水 量 实测 
值 相关 性 分 析 
Fig. 5 Correlation analysis between volumetric soil water 

content predicted by 4 models with5 depth ranges and drying 
method observations with 0 to 50 cm under arid soil and moist 


soil conditions 
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深度 的 增加 而 逐渐 增 大 ， 相 关系 数 在 0~50 cm 深 
度 范围 达到 最 大 值 (R=0.626); 其 他 三 个 模型 
的 相关 性 变化 趋势 相同 ， 均 在 0~40 cm 深度 范围 
达到 最 大 值 ， 其 中 了 Roth 模型 相关 性 最 好 (R= 
0.645)。 而 4 个 模型 在 0~10 cm Ail 0—20 cm 深度 
范围 的 决定 系数 较 低 (R 在 0.3 上 下 浮动 ), 在 
0—30 cm 深度 范围 开始 显著 增加 。Roth 模型 在 
0 一 10 cm 至 0—40 cm 深度 范围 均 具有 最 小 的 
RMSE, mÆ 0—50 cm 深度 范围 内 Herkelrath 模 
型 RMSE 最小。 在 土壤 湿润 条 件 下 ，4 个 模型 相 
关 性 变化 趋势 相同 Roth 模型 在 O~40 cm 深度 
范围 相关 性 最 好 (R= 0.750) ， 并 且 在 0—10 cm 
28 0—30 cmRMSE 最 小 。 

在 有 效 反 演 深度 范围 内 ,通过 对 比 两 种 土壤 
水 分 条 件 下 4 个 模型 反 演 值 与 实测 值 的 RMSE 平 


0—30 cm 深度 范围 AE 最 小 ， 而 在 O~40 cm 和 
0—50 cm TEE T [S], Roth 模型 的 反 演 值 与 实测 值 
相 比 偏 高 。 

此 外 ， 随 着 土壤 采样 深度 的 逐渐 增加 ，RE 
逐渐 降低 ( 表 3)。 在 无 降雨 条 件 下 ，Roth 模 型 
在 0~10 cm 至 0~40 cm 范围 内 RE 最 小 ; 在 0~ 
50 cm REWENA, Topp HERE WHK, HES 
个 模型 RE 值 较 接近 。 在 土壤 湿润 条 件 下 ，4 个 
模型 的 RE 均 在 0~40 cm 深度 范围 内 达到 最 小 值 ; 
Roth 模型 在 0—10 cm 至 0~30 cm 深度 范围 内 与 
其 它 三 个 模型 相 比 RE 最 小 ， 而 在 0~40 cm 至 0 一 
50 cm 深度 范围 内 ， 其 RE 较 大 。 

表 3 模型 反 演 值 与 土壤 含水 量 实测 值 的 相对 误差 (n=30) 


Table 3 The relative error between model inversion value 


and volumetric soil water content (n=30) 


均值 ( 表 2)，Roth 模 型 在 5 个 深度 范围 内 的 — ER 
RMSE 平均 值 明 低 于 其 它 3 个 模型 . 表 明 Roth Topp Roth Herkelrath Ferre 
0—10 37:3 26.3 35.9 36.7 
SE p Jat YR HE t FE H < Ju 
模型 在 有 效 反 演 深度 范围 内 具有 更 稳定 的 注 om so my 2 xs 
a LE 
REJJ o TART 030 163 107 172 167 
7 旱 条 件 
表 2 不 同 深度 范围 模型 反 演 值 与 土壤 含水 量 实测 值 的 均 4 40 5.8 3.6 7.6 65 
方 根 误差 平均 值 (n=30) 0~50 3.8 3.0 2.8 3.4 
Table 2 The mean RMSE between model inversion value 0—10 17.8 15.1 19.3 184 
and volumetric soil water content in different depth ranges 土壤 偏 湿 人 ae ee ane em 
B | 0-30 10.3 5.9 11.6 10.9 
(n=30) 润 条 件 
0 一 40 3 4.6 4.3 3.9 
均 方 根 误差 平均 值 RMSE/(m m?) 0~50 5.3 5.4 5.1 5.1 
土壤 水 分 条 件 
Topp Roth Herkelrath Ferre Bk I "T 
te Zh M 
Te 0.0563 — 0.0401 0.0559 0.0561 因此 ， 基 于 在 两 种 土壤 水 分 条 件 下 的 41 
i FE Ee y 与 实测 A X VEZ 
oo 00335 Wo oo 模型 反 演 值 与 实测 值 的 相关 性 分 析 及 包括 


图 6 是 两 种 土壤 水 分 条 件 下 ， 有 效 反 演 深度 
范围 内 模型 反 演 值 与 实测 值 的 AE 变化 趋势 图 。 
从 反 演 结果 来 看 ，4 个 模型 的 反 演 值 普遍 低 于 实 
测 值 。 在 无 降雨 条 件 下 ， 随 着 深度 的 增加 ， 模 型 
反 演 误差 逐渐 降低 。 在 0~10 cm 至 0~40 cm T 
度 范围 内 ， 与 其 它 模 型 相 比 Roth 模 型 反 演 值 的 
AE 较 低 。0~-50 cm 深度 范围 内 ，4 个 模型 的 AE 
较 接近 ，Herkelrath 模 型 AE 最 低 ，0.02 m/m, 
在 地 表 湿 润 条 件 下 ，4 个 模型 均 在 0~40 cm 深 
度 范围 误差 最 小 。 其 中 Roth 模 型 在 0~10 cm 至 


RMSE 在 内 的 3 种 误差 分 析 结 果 ， 本 研究 选 定 
1000 MHz 探 地 雷达 在 试验 田 块 的 最 优 模 型 为 
Roth 模 型 。 在 无 降雨 条 件 下 ， 最 佳 反 演 深度 为 
50cm; 在 土壤 湿润 条 件 下 ， 最 佳 反 演 深度 为 
40 cm. 


3.3 相 邻 反射 层 间 土壤 含水 量 反 演 精 度 
利用 雷达 波 通 过 不 同 深度 反射 层 的 速度 值 和 
Dix 公式， 获取 雷达 波 在 10~20 cm, 20~30 cm, 
30—40 cm Fil 40~50 cm 土壤 层 的 层 速度 ， 并 利用 
公式 (2) 获取 不 同 层 间 土壤 的 相对 介 电 常数 ， 
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d 
o 


0.0 


测 值 的 绝对 误差 (mm ) 


试 -10|. 


1 
2 


I 
a 


模型 反 演 值 与 实测 值 的 绝对 误差 /mzm3) 


模型 反 演 值 与 


0~10 0-20 0-30 0-40 0-50 


0-10 0-20 0-30 0-40 0-50 
深度 /cm 深 


TDR 在 实验 室 中 建立 ， 对 于 本 次 试 
验 结果 造成 较 大 误差 可 能 是 其 对 于 
农田 环境 的 层 间 土 壤 水 分 反 演 的 适 
用 性 不 强 所 致 。 


4 讨论 


本 研究 通过 设置 探 地 雷达 CMP 
法 田间 土壤 水 分 监测 试验 ， 期 望 控 
究 探 地 雷达 对 不 同 深度 及 相 邻 反映 
层 间 土 壤 仿 水量 的 反 演 精度 。 分 析 


度 / ———nÜ iei 
zm 雷达 波 速 度 及 介 电 常数 ,结合 


(a) 无 降雨 土壤 偏 干旱 条 件 土 壤 含 水 量 (b) 降 雨 后 土 壤 偏 湿润 条 件 土 壤 含 水 量 模 RMSE R AE All RE 确定 田 块 尺 


模型 反 演 值 与 实测 值 绝对 误差 分 析 


Fig.6 Absolute error between volumetric soil water content ( VSWC) predicted 


by 4 models with 5 depth ranges and drying method observations with 0 —50 cm 


under arid soil and moist soil conditions 


再 结合 Roth 模 型 获取 不 同 深度 相 邻 反射 层 间 土 
壤 售 水 量 的 模型 反 演 值 。 模 型 反 演 值 与 实测 值 的 
精度 验证 如 表 4 所 示 ， 无 降雨 条 件 下 20~30 cm 
深度 范围 误差 最 小 ， 其 余 深 上 度 范围 误差 较 大 ， 
RE fE 40—50 cm 深度 范围 最 高 可 达 41.3%。 降 
雨 后 土 壤 湿 润 条 件 下 ，4 个 深度 范围 中 仅 在 
40—50 cm 反 演 精度 较 高 ， 并 且 10~20 cm £ 30— 
40 cm 深度 范围 的 反 演 精度 随 着 深度 的 增加 ， 误 


型 反 演 值 与 实测 值 绝 对 误差 分 析 
图 6 无 降雨 偏 干旱 和 降雨 后 偏 湿 润 土壤 条 件 下 0—50 cm 5 个 深度 范围 4 种 
模型 预测 值 与 烘 干 法 土壤 体积 含水 量 实测 值 的 绝对 误差 


度 最 优 模 型 、 雷 达 波 有 效 反 演 深度 
以 及 两 种 土壤 水 分 条 件 下 的 最 佳 反 
演 深 度 。 

分 析 雷 达 波 对 不 同 深度 土壤 水 
分 的 监测 精度 ， 反 射 波 对 土壤 含水 
量 的 反 演 精度 高 于 地 表 波 。 可 能 
表 4 两 种 土壤 水 分 条 件 下 0~50 om 深度 范围 相 邻 反射 

层 间 土壤 含水 量 反 演 结果 (n=30) 


Table 4 Inversion results of soil water content between adja- 


cent reflective layers in 0—50 cm depth range under two 


soil water conditions (n=30) 


HBAR, HBAS O~10 cm 至 0—40 cm 深 
度 范 围 的 反 演 结果 相反 。 

因此 ， 依据 本 次 试验 结果 可 以 看 出 ， 利 用 探 
地 雷达 及 Dix 公式 对 不 同 深度 相 邻 反射 层 间 土壤 
水 分 的 反 演 具 有 局 限 性 ， 其 反 演 精度 有 待 进一步 


深度 范围 / OR gc 3) RE/% Bo 
ilt5icm msm 
^ 条 件 ^o fmm?) 
MFE -0.38 12.7 0.0402 
10~20 c 
偏 湿润 -0.67 18.9 0.0563 
MTE 0.16 9.3 0.0308 
20~30 
ETAR 0.68 28.3 0.0806 
MFE 0.78 19.2 0.0611 
30~40 uu 
偏 湿润 1.40 43.7 0.1332 
偏 干旱 1.30 41.3 0.1347 
40~50 TM 
DETAR 0.15 19.5 0.0126 


加 强 。 无 降雨 条 件 下 ， 对 层 间 土壤 的 反 演 精度 更 
高 ;然而 ， 对 于 从 地 表 至 反射 层 间 土壤 水 分 的 反 
演 精度 来 说 ， 无 降雨 条 件 的 反 演 精度 低 于 降雨 后 
土壤 湿润 条 件 下 的 反 演 精度 。 由 此 可 见 ， 降 雨 对 
探 地 雷达 的 反 演 精度 产生 一 定 的 影响 ， 其 原因 一 
方面 可 能 是 由 于 降雨 后 表层 土壤 水 分 未 完全 下 渗 
造成 反 演 误差 ; 9 — M, Roth 模型 是 基于 


于 试验 田 块 表 层 土壤 受 耕 作 过 程 等 影响 ， 其 土壤 
状态 易 被 破坏 ， 导 致 地 表 0~10 cm 土壤 样本 实测 
值 与 地 表 波 反 演 值 出 现 较 大 偏差 。 通 过 观察 两 种 
土壤 水 分 条 件 下 的 雷达 图 像 ， 试 验 田 块 深层 土壤 
有 明显 的 反射 层 ， 并 且 空 间 变 异性 较 小 ， 表 明 双 
曲线 拟 合法 可 以 得 出 较 准 确 的 波 速 值 。 但 是 雷达 
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图 像 能 否 存 在 明显 的 反射 波 取决 于 地 理 位 置 及 土 
壤 类 型 ， 前 人 人 研究 表明 在 没有 明确 地 下 反射 层 的 
条 件 下 ， 地 表 波 法 仍然 是 测定 浅 表层 土壤 含水 量 
的 有 效 方法 ,今后 应 加 强 探 地 雷达 在 不 同 土 
壤 类 型 的 有 效 反 演 深度 探究 和 地 表 波 的 反 演 
精度 。 

土壤 含水 量 模型 精度 影响 探 地 雷达 的 反 演 精 
度 ， 本 次 试验 得 出 Roth 模 型 是 田 块 尺度 最 优 反 
演 模型 。 但 Roth 模 型 的 参数 是 基于 时 域 反 射 法 
和 土壤 样本 于 实验 室 建立 ,与 农田 环境 差别 较 
大 。 在 不 同 试验 区 域 若 想 提 高 探 地 雷达 的 反 演 精 
度 ， 需 要 依据 试验 区 实际 土壤 水 分 情况 对 Roth 
模型 参数 进行 校正 。 


5 结论 


为 确定 田 块 尺度 下 探 地 雷达 对 不 同 深度 及 相 
邻 反 射 层 间 土壤 含水 量 的 反 演 精度 、 最 优 反 演 模 
型 、 有 效 反 演 深度 及 最 佳 反 演 深度 ， 本 研究 在 无 
降雨 偏 干旱 和 降雨 后 偏 湿润 两 种 土壤 水 分 条 件 下 
各 选取 15 个 测 点 , 采用 1000 MHz 中 心 频率 探 地 
雷达 ， 进 行 CMP 法 测定 土壤 体积 含水 量 。 将 4 种 
模型 反 演 值 与 烘 干 法 实测 值 进行 精度 验证 后 得 出 
结论 : 

(1) 1000 MHz 探 地 雷达 的 有 效 反 演 深度 范 
围 为 0~50 cm; 在 土壤 偏 干旱 条 件 下 最 佳 反 演 深 
度 为 50 cm， 在 土壤 湿润 条 件 下 最 佳 反 演 深度 为 
40 cm. 

(2) 与 其 它 3 个 模型 相 比 ，Roth 模 型 相关 性 
Hf. Rei 0.750; 且 反 演 土 壤 含水 量 值 更 
稳定 ， 其 中 在 土壤 偏 干 旱 和 偏 湿润 条 件 下 RMSE 
平均 值 分 别 为 0.0401 和 0.0335 m/m; RE 最 低 为 
3.0%. 

(3) 利用 探 地 雷达 及 Dix 公式 对 不 同 深度 相 
邻 反 射 层 间 土壤 水 分 的 反 演 具有 局 限 性 ， 其 反 演 
精度 有 待 进一步 加 强 。 
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Monitoring Specified Depth Soil Moisture in Field Scale 
with Ground Penetrating Radar 


ZHANG Wenhan', DU Keming", SUN Yankun", LIU Buchun’, SUN Zhongfu’, 
MA Juncheng’, ZHENG Feixiang* 


(1. College of Resources and Environment, Northeast Agricultural University, Harbin 150030; 2. Institute of Envi- 
ronment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081) 


Abstract: Ground-penetrating radar (GPR) is one of the emerging technologies for soil moisture measurement. However, the 
measurement accuracy is difficult to determine due to some influence factors including radar wave frequency, soil texture type, 
etc. The GPR equipment with 1000 MHz center frequency and the measurement method of common midpoint (CMP) were ad- 
opted in the research to collect radar wave raw data in the selected field area under arid soil and moist soil conditions. The trans- 
mitter and receiver antennas of the GPR equipment were moved 0.01 m respectively in opposite directions on each radar wave 
raw data collection. Therefore, a CMP radar image consisted of 100 pieces of radar wave raw data by increasing the antenna dis- 
tance from 0 m to 2 m. Each radar wave raw data indicated that the radar waves were reflected in the reflective layer with differ- 
ent dielectric constant under the same antenna distance. And the reflected and refracted radar waves were acquired by the receiv- 
ing antenna at different two-way travel time respectively, and recorded in the computer. The collection of CMP soundings aimed 
to determine the inversion accuracy, optimum inversion depth, effective inversion depth and optimal inversion model of soil 
moisture content at different depth ranges and adjacent reflective layers by GPR at field scale. The reflected and refracted radar 
wave data were extracted from the raw data. The velocities of the surface waves and reflected waves were obtained respectively 
from the line slope of the surface wave data and the hyperbolic curves fitting of the reflected wave data. In addition, the relative 
dielectric constant of the soil at specified depth were deduced according to the soil dielectric constant and its reflected wave ve- 
locity. Moreover, 4 different models including Topp, Roth, Herkelrath and Ferre were used to figure out the soil volumetric wa- 


ter content inversion. Meanwhile, the measured data of soil volumetric moisture content obtained by oven drying method were 


used to verify the accuracy of the inversion results. The results showed that the effective inversion depth of 1000 MHz GPR 


ranged from 0 to 50 cm. The best inversion depth was 50 cm in arid soil and 40 cm in moist soil. The Roth model had the best 
correlation and stability with the highest R’ was 0.750, the Root Mean Square Error (RMSE) was 0.0114 m’/m’ and the lowest 
Relative Error (RE) was 3.0%. The GPR could possess the capacity of quick, precise and non-destructive measurement of speci- 
fied depth soil moisture in field scale. The inversion model of soil moisture content needs to be calibrated according to different 
soil conditions. 

Key words: ground penetrating radar; soil water content monitoring; common midpoint; optimum inversion depth; optimum in- 


version model; Roth model 
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